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Results  on  the  detection  of  narrow-band 
low-level  signals  In  a  ship-induced  noise 
field  using  high-resolution  beainformiog 

S.  Jesus  and  R.M.  Heitmeyer 

Abstract:  In  tnis  report  the  performance  of  a  high- resolution  beamformer 
is  compared  to  that  of  the  conventional  beamformer  in  detecting  a  low-level 
narrow-band  signal  in  the  vicinity  of  a  number  of  high-level  noise  sources 
This.  scenario  may  arise  in  areas  whete  shipping  noise  is  predominant,  the 
acoustic  held  being  dominated  to  a  large  part  by  the  contributions  of  in¬ 
dividual  ships  In  this  study  a  comparative  detection  performance  test  is 
conducted  with  <1  simulated  source  embedded  in  a  real  ship-induced  noise 
held  and  received  on  a  64-element  towed  array  The  results  shown  sug¬ 
gest  that  in  such  noise  fields  the  detection  performance  of  the  conventional 
pianewave  beamformer  is  greatly  reduced  due  to  a  lack  of  resolving  power 
High-resolution  techniques  may  overcome  this  difficulty  and  a  detection  im¬ 
provement  of  5%  to  20%  is  obtained  in  the  detection  of  a  low-level  target  in 
the  frequency  range  below  250  Hs.  When  the  signal-to-noise  ratio  increases 
the  performance  of  both  algorithms  is  comparable. 
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1.  Introduction 

For  a  conventional  planewave  beamforilier  an  increase  in  the  array  length  results  in 
an  increase  in  the  signal-to-noise  ratio  gain  and  hence  an  increase  in  the  detection 
performance  If  the  noise  field  is  formed  by  a  single  point  source  in  spatially  white 
gaussian  noise  then  the  conventional  beamformer  is  the  optimum  beamfonner  and 
the  increase  in  the  signal-to-noise  ratio  gain  is  the  only  benefit  of  increased  aperture 
length 

For  ship-induced  noise  fields,  where  the  noise  is  largely  due  to  the  contribution  of  in¬ 
dividual  ships,  an  increase  in  the  array  length  also  results  in  an  increase  in  the  extent 
to  which  the  contributions  from  the  individual  ships  are  resolved  in  steering  angle 
Thus,  for  ship-induced  noise  fields,  increasing  the  aperture  length  affects  detection 
performance,  not  only  by  increasing  the  array  gain  but  also  by  increasing  the  extent 
to  which  noise  sources  are  resolved.  In  such  noise  fields  the  limited  resolution  power 
of  the  conventional  beamformer  may  lead  to  a  situation  where  low-level  sources  may 
be  masked  by  nearby  high-level  sources.  Over  the  past  two  decades  a  number  of 
high-resolution  spatial  array  processing  algorithms  have  been  developed  to  improve 
the  ability  of  sonar  systems  to  angularly  resolve  closely-spaced  sources  [1-5]  For 
ship-induced  noise  fields  these  algorithms  may  offer  an  alternative  to  increasing  the 
aperture  length  An  overview  and  comparative  simulated  study  of  the  best-known 
high-resolution  techniques  applied  to  underwater  acoustics  can  be  found,  for  exam¬ 
ple,  in  [6],  The  Ambient  Noise  Group  (ANG)  at  SACLANTCEN  has  been  involved 
in  developing  and  testing  such  techniques  to  improve  angular  resolution  of  short 
linear  towed  arrays  for  horizontal  noise  directionality  measurements  [7,8] 

In  this  report,  we  present  an  experimental  evaluation  of  the  detection  performance 
of  a  high-resolution  beamformer  compared  to  that  of  the  conventional  planewave 
beamformer.  The  high-resolution  technique  used  in  the  evaluation  is  the  WB5  algo¬ 
rithm  developed  by  J.-L.  Berrou  as  a  modification  of  the  WB*  described  in  Wagstaff 
et  al.  [8].  Our  choice  of  the  WB*  algorithm  was  largely  motivated  by  data-processing 
considerations  rather  than  performance  considerations.  The  results  are  described  by 
receiver  operating  characteristics  (ROC)  curves  (probability  of  detection  vs  proba¬ 
bility  of  false  alarm)  obtained  both  for  the  high-resolution  and  conventional  beam- 
formers. 

The  data  used  for  the  experimental  detection  evaluation  is  the  high-resolution  data¬ 
base  acquired  during  the  ANG  cruise  1/84  in  the  Gulf  of  Genova.  The  choice  of  this 
working  area  was  motivated  by  the  high  concentration  of  merchant  ship  traffic.  This 
database  provides  a  relatively  long  time  record  of  data  from  a  64-hydrophone  towed 
array  with  a  known  controlled  multitone  source  signal  embedded  in  ship- dominated 
noise. 
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2.  Description  of  the  algorithms 


The  detection  performance  evaluation  is  obtained  by  means  of  a  conventional  detec¬ 
tion  algorithm  operating  on  a  time  series  of  frequency-wavenumber  estimates  com¬ 
puted  by  either  the  conventional  beamformer  or  the  high-iesolution  beainformer 
Each  frequency-  wavenumber  estimate  is  in  turn  computed  in  the  conventional  man¬ 
ner  as  a  time-average  of  the  two-dimensional  FFT  of  the  array  sensor  outputs 
Thus  the  M-averaged  frequency  -wavenumber  estiniate  at  discrete  frequency  /*  and 
wavenumber  at  time  t  is  computed  according  to 

N 

B(A,Kn.t)  =  M  ~ 1  (la) 

where  the  ‘unavt.aged’  frequency-wavenumber  estimate  at  time  tm  is  given  by 


-  1  J 

52  MW(/*.*m)exp(;2ir*n.') 
i=s 


(lb) 


where  yi(/k,tm)  is  the  complex  amplitude  of  the  1th  hydrophone  output  at  fre¬ 
quency  fk  and  time  'snapshot'  tm,  h(l)  is  the  spatial  shading  function,  and  Mb  is 
the  number  of  beams  The  complex  hydrophone  amplitudes  are  obtained  from  a 
time-to-frequency  fourier  transform  on  the  data  acquired  over  the  recording  interval 
[rm  -  Tr,tm].  The  discrete  frequencies  /*  with  fc  =  1, .  .  . ,  X  take  on  equally-spaced 
values  in  the  frequency  band  of  interest  with  a  frequency  spacing  A  /  =  Tt~ 1 .  The 
discrete  normalized  wavenumbers  k„  take  on  the  values 

*„  =  n/Mb,  n  -  1, . . . ,  Mb  -  l,  (2) 

where  the  number  of  beams  Mb  is  an  integer  multiple  of  the  number  of  hydrophones 
A'p.  If  M b  >  M p  the  spatial  shading  function  is  h(l)  =  0  with  l  =  Nv  +  1,.  ..Mb 
(zero  padding) 
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3.1  CONVENTIONAL  BEAMEORMER 


The  conventional  beamformer  operates  on  the  array  by  coherent  summation  of  the 
N p  sensor  outputs  (for  a  given  frequency  fu)  with  the  appropriate  phase  shifts 
to  steer  a  'beam’  n  to  direction  9„.  This  is  done  by  identifying  the  normalieed 
wavenumbers  of  (2)  with  the  beam  steering  angles  defined  by  the  relationship 


*n 


=  arcsin 


(3) 


where  9„  €  [-90°, 90°)  and  k),  =  (n  -  6  [0,1)  for  n  =  1,  ..  ,/Vb  if  the  array 

design  frequency  /d  -  c/2d  >  /*  which  is  the  case  in  this  report. 

•  For  fu  -  It  the  direction  9t  =  -90°  corresponds  to  the  aft  endfire  beam,  the 
direction  0(jvb  +  i)/J  =  arcsin( )  corresponds  to  the  broadside  beam  and 
the  direction  -  arcsin(l  -  2A'b‘’1 )  corresponds  to  the  near  forward  endfire 
beam. 


•  For  fu  <  ji  Ec.  (3)  gives  sin9„  outside  of  [-1, 1],  the  corresponding  beam 
wavenumbers  are  not  related  to  a  real  acoustic  direction  and  are  commonly 
termed  virtual  beams;  conversely  those  lying  in  the  real  acoustic  region  are 
termed  real  beams 


3.3.  HIGH  RESOLUTION  BEAMFORMER 

In  most  high-resolution  methods  some  a  pnort  information  is  required  on  the  form 
of  the  model  for  the  field  generating  process  and/or  for  the  estimation  of  the  model 
parameters.  Mismodeling  and  sub-optimum  parameter  estimation  are  very  often 
confused  in  the  real  case 

To  avoid  these  problems  a  different  solution  was  proposed  by  Wagstaff  et  a).,  the 
Wagstaff-Berrou  broadband  (WB1)  algorithm  (8).  In  this  report  we  shall  use  a 
modified  version  of  this  algorithm  due  to  J.-L  Berrou  (WB*);  for  convenience  the 
algorithm  equation  flowchart  is  given  in  the  appendix.  Here  we  will  simply  'ecall 
the  basic  ideas  of  the  WB  techniques 

In  the  WB  algorithms  noise  field  estimetion  is  considered  as  an  inverse  problem 
with  the  conventional  beamforraer  estimate  as  input.  The  WB  algorithms  are  based 
on  the  assumption  that  the  acoustic  field  incident  on  the  array  can  be  described 
as  a  sum  of  uncorrelated  plane  waves  with  wavenumber  values  sc,,  i  —  1,.  A^,. 

Thus  the  mean  of  the  conventional  beamformer  output  at  time  t,  frequency  fu  and 
wavenumber  sc„  can  be  expressed  by  the  discrete  wavenumber  convolution  relation¬ 
ship 

Nk-I 

B(/k,n.f)=£(B(/*,«„,t)]  =  AT-'  Y.  »  =  1 . Nb,  (4) 

«  =  0 


,  w 

I 
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where  F<  -  F(f*.K.i,t)  with  t  =  l,  . ,  .Vh  is  the  discrete  wavenumber  spectrum 
given  by  a  hypothetical  array  of  A'h  >  Nb  sensors  with  the  same  hydrophone  spacing 
as  the  actual  array  and  Pn.,  is  the  wavenumber  spectrum  of  the  spatial  shading 
function  h(l),  i.e 


."a-* 


F,  -  ^  M0«pI>2»k.1| 


(5) 


evaluated  at  the  point  n  i  The  basic  idea  (knowing  the  array  beanipattern  Pi) 
is  to  deconvolve  the  discrete  wavenumber  spectrum  F,  from  the  measured  conven¬ 
tional  beamformer  bearnnumber  spectrum  B(/k,n,t).  In  the  WB  algorithms  all  the 
pectrum  quantities  are  expressed  in  a  decibel  scale  which  emphasizes  the  impor¬ 
tance  of  the  low  acoustic  levels  and  facilitates  the  estimation  of  low-level  sources  in 
the  presence  of  high-level  noise  sources  (8|.  The  computer  implementation  of  the 
algorithm  is  somewhat  simplified  by  this  decibel  notation  but  the  inverse  problem 
turns  out  to  be  highly  nonlinear  Also,  as  the  number  jVh  of  discrete  values  of  the 
field  F,  it  much  larger  than  the  number  of  beams  the  inversion  of  (4)  leads  to 
an  indetermination  of  minimum  order  -  Nt,. 


The  WB  algorithms  give  an  approximate  solution  to  this  problem  by  means  of 
a  gradient  iterative  technique  based  on  the  minimization,  at  each  iteration  k,  of 
the  mean  square  eror  between  the  predicted  beamformer  response  (4)  for  the  Jtth 
iteration  discrete  wavenumber  spectrum  F,1  Rnd  the  actual  measured  beamformer 
response. 

The  main  differences  between  the  original  algorithm  ( WB1 )  and  the  actual  algorithm 
(WB3)  consists  in  the  first  guess  of  the  field  Ff  ami  the  correction  applied  at  each 
iteration  to  pass  from  F,*1  to  F1I,+  I  The  iteration  process  is  stopped  either  when  the 
mean  rquare  error  passes  below  a  given  threshold  or  when  a  given  fixed  number  of 
iterations  has  been  reached 


The  WB1  high  resolution  algorithm  was  qualitatively  shown  to  give  similar  results  to 
the  majority  of  well-known  high- resolution  techniques  [8].  A  simulated  comparative 
quantitative  study  of  resolution,  detection  and  accuracy  performance  is  given  in  [9). 
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3.  Methodology  of  detertion 

The  purpose  of  this  section  is  to  describe  the  detection  procedure  adopted  to  carry 
out  a  statistically  significant  signal  detection  study  comparing  the  conventional 
planewave  heamformer  and  the  high-resolution  beamformer  on  a  real  ship-induced 
environment. 

This  study  is  restricted  to  the  test  of  presence  or  absence  of  a  single  najrow  band 
signal  (target)  when  the  direction  and  frequency  signature  are  known.  Standard 
detection  procedures  normally  use  a  frequency-asiniuth  display  and  the  assumption 
that  the  frequency  bins  apart  from  the  signal  frequency  are  noise-only  frequencies. 
These  noise-only  frequency  bins  are  then  averaged  to  produce  an  estimate  of  the 
mean  noise  power  In  this  case,  the  detection  problem  leads  to  a  simple  conventional 
testing  against  a  given  power  threshold  based  on  the  previous  estimated  mean  noise 
power 

For  the  database  discussed  in  this  report,  frequency-wavenumber  estimates  are  ob¬ 
tained  for  266  discrete  frequencies  ranging  from  112  5  Hi  to  760  Hs  at  2  5  Hi  incre¬ 
ments  Each  surface  is  obtained  for  an  averaging  time  of  slightly  less  than  3  min; 
Nb  =  128  beams  were  formed  which  is  twice  the  number  of  hydrophone  elements; 
H&nn  shading  was  used  both  in  time-to-frequency  and  hydrophone-to-wavenumber 
transforms.  For  testing  purposes  a  multitone  controlled  source  emitted  from  a  known 
location  about  75  n.mi  away  from  the  measurement  site  was  used  as  a  target.  The 
controlled  source  input  spectrum  is  formed  by  a  number  of  CW  lines  ranging  from 
133.33  Hs  to  600  Hs  with  a  frequency  increment  of  33.33  Hi  and  variable  power  level 
A  reference  peak  is  located  at  the  frequency  of  500  Hr  which  is  used  as  reference  for 
source  asimuth  estimation.  The  five  lowest  frequency  lines  were  used  as  test  signals 
in  this  study  (see  Table  1). 

Let  us  define  as  'signal  frequencies'  the  discrete  frequencies  of  the  frequency -wave- 
number  spectrum  estimates  that  are  close  to  the  source  CW  line  frequencies  (see 
Table  1).  By  assumption  the  other  frequencies  will  be  called  the  'noise-only  frequen 
cies’. 

For  detection  and  false  alarm  probabilities  estimation  a  detection  is  obteuned  when¬ 
ever  a  line  is  detected  at  a  ‘signal  frequency’,  and  a  false  alarm  is  obtained  whenever 
a  line  is  detected  at  a  'noise-only  frequency'  Thus,  the  empirical  estimates  of  the 
detection  probability  Va  and  the  false  alarm  probability  Vt ,  are  respectively  given 
by 


Pa 


number  of  line  detections  at  signal  frequency 
total  number  of  estimates 


(6a) 
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Table  1 

Transmitted  sourer  signal  cliatacieiislica 


Source 

frequencies  (H») 

Relative 
power  (dB) 

Discrete 

frequencies  (Hi) 

133.33 

55.28 

132  50 

166  66 

40.23 

167.50 

200  00 

42  60 

200  00 

233.33 

-40  30 

232.50 

266  68 

-36  10 

287.50 

~  number  of  line  detections  at  noise-only  frequency 

>7  ,  -  - - - - - r* - — - : —  ■  fob) 

total  number  of  estimates 

For  the  conventional  planewave  bearnformer  estimate  one  can  define  the  simplest 
detection  test  by 

(7> 

where  fk  may  be  a  signal  or  a  noise-only  frequency,  depending  upon  whether  we 
estimate  the  detection  probability  or  the  false  alarm  probability;  6>,  is  the  source 
asimuth  beam  steering  angle  and  lj  is  the  discrete  time  at  which  the  estimate  was 
obtained 

For  the  high-resolution  WB1  bearnformer  the  detection  test  is  defined  by  the  expres¬ 
sion 

MF(/*,0.,f,)>  (8) 

where  MF(/i,,  t;)  is  defined  as  the  maximum  level  (expressed  in  dB)  of  the  WB3  field 
estimate  within  the  width  of  the  conventional  planewave  bearnformer  mainlobe  at 
the  considered  discrete  frequency  fk  and  for  the  source  signal  beam  steering  angle 
6,.  The  detection  threshold  is  simply  7 ,  expressed  in  dB. 

In  (7)  7j  is  the  mean  noise  power  ‘background’  dependent  threshold  at  time  t}.  In 
practice  the  noise  power  ‘background’  may  be  estimated  either  by  an  average  of  a 
number  of  ‘noise-only  frequency’  bins  in  the  signal  direction  or  by  an  average  of  the 
noise  levels  in  contiguous  beam  steering  directions  at  a  noise-only  frequency,  or  even 
some  combination  of  both  estimates.  In  the  general  case  the  detection  threshold  is 
given  by 

1  1  Q 

7j  =  7T  +  7 (9) 
v  1=1  *=! 

where  L  is  the  number  of  noise-only  frequency  bins  and  Q  is  the  number  of  contiguous 
beams;  77  is  a  variable  detection  threshold  depending  on  the  false  alarm  rate  for  a 
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Used  signal- to- noise  ratio.  In  (9)  the  symbol  C  stands  for  the  conventional  or  the 
WB*  high- resolution  beomformer  level  (at  the  specified  frequency,  steering  angle 
and  time)  depending  on  whether  a  conventional  or  high-resolution  detection  test  is 
being  performed. 

The  results  are  shown  in  the  form  of  a  series  of  receiver  operating  characteristics 
(ROC)  curves,  i.e.  estimated  Vi  as  a  function  of  estimated  Vt  4  for  variable  detection 
threshold  77,  for  the  lowest  signal  frequencies  of  the  source  signed  (see  Table  1). 
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4.  Results 

The  high-resolution  database  used  in  this  study  is  composed  of  about  6  h  of  array 
data  recording  which  provides  a  total  of  120  frequency-wavenumber  estimates.  The 
sampling  frequency  is  F,  =  2560  Hs  for  an  FFT  block  size  of  .V,  -  1024.  The  number 
of  averages  (time-bandwidth  product)  is  ,V  60  Figure  1  shows  an  example  of  a 
conventional  beamfornier  output  (a)  and  the  respective  WB  high-resolution  estimate 
(b)  for  a  frequency  /,  -  200  Hz. 


Fig.  1.  Conventional  beamfoimer  (a)  and  WB'  high-resolution 
(b)  estimates;  number  of  beams  .Vb  =  129  (a)  and  iVb  ^  512  (b); 
number  of  averages  X  —  150,  frequent)  /,  ^  200  Hz 


Results  are  presented  as  a  series  of  five  empirical  ROC  curves  (percentage  of  detec¬ 
tion  ns  percentage  of  false  alarms)  for  a  threshold  variation  ymjn  <  7T  <  7n,.b 

7„.,„  =  0.316  (  -  -5  dB), 

7m.,  -2  81  (=  4.5  dB) 

There  are  a  total  number  of  40  discrete  points  on  each  ROC  curve.  In  practice  the 
values  L  ~  3  and  Q  -  3  were  chosen  for  the  expression  of  the  mean  noise  power 
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estimator  of  (10).  The  three  noise-only  frequencies  were  chosen  one  below  the  signal 
frequency,  one  between  the  signal  and  the  noise  frequencies,  and  one  above  the  noise 
frequency.  The  five  signal  frequencies  in  Table  1  (132.5,  167  5,  200.0,  232.5  and 
267.5  Hr)  have  been  processed  giving  the  results  shown  in  Figs.  2-6.  To  facilitate 
comparison  the  ROC  curves  of  both  the  conventional  planewave  beamformer  and 
the  high-resolution  beamforme,-  are  plotted  on  each  figure. 

For  the  four  lowest  frequency  lines  (/,  =  132.5,  167.5,  200  and  232.5  Hr)  the  con¬ 
ventional  beamformer  showed  a  relatively  poor  performance,  only  slightly  better 
than  the  chance  diagonal  defined  by  Pd  =  F\  For  these  frequencies  the  detection 
performance  of  the  high-resolution  beamformer  is  from  5%  to  20%  larger  than  that 
achieved  by  the  conventional  beamformer  Note  that  this  improvement  was  obtained 
for  the  range  of  false  alarm  rates  of  most  interest  to  us,  i.e.  for  5%<  P(  <40%. 

For  the  highest  level  line  (/,  —  267  5  Hz)  the  performance  of  the  two  beamformers  is 
comparable,  attaining  relatively  high  detection  probability  rates  for  the  same  false 
alarm  probability  range. 

These  results  suggest  that  the  poor  performance  of  the  conventional  plane-wave 
beamformer  for  the  lowest  frequencies  is  due  not  only  to  the  low  signal-to-noise 
ratio  of  the-  signal  under  detection  but  also  to  a  lack  of  resolving  power.  This  lack 
of  resolving  power  could  be,  at  least  partially,  recovered  by  using  a  high-resolution 
beamformer  giving  rise  to  an  increase  in  detection  performance  without  increasing 
the  array  aperture. 

When  operating  at  the  highest  frequencies,  where  the  signal  to  noise  ratio  and  the 
resolving  power  are  higher,  the  two  beamformers  gave  the  same  detection  perfor¬ 
mance. 


SaCLANTCEN  SR-144 


Fig.  2.  Empirical  ROC  curves  of  the  conventional  and  high- 
resolution  beamformer  for  frequencies  /,  =  132.5  H*  and  /B  = 
150  He 


Fig  3.  Empirical  ROC  curves  of  the  conventional  and  high- 
resolution  beamformer  for  frequencies  /,  =  167.5  H»  and  /«  = 
182  5  He 
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Fig  6  Empirical  ROC  curvee  of  the  conventional 
resolution  beamformer  (or  frequencies  /,  =  267.5  Hs 
282.5  Hs. 
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S.  Conclusion 

In  this  report  the  results  of  the  detection  of  &  low-level  narrow-band  signal  in  real 
coherent  ship-induced  noise  are  presented.  The  detection  performance  behaviour 
of  both  the  conventional  beamformer  and  a  high-resolution  beamformer  (WB 
gorithm)  are  shown  under  the  form  of  ROC  curves  for  different  frequencies  and 
signal-to-noise  ratio. 

The  results  obtained  on  this  particular  set  of  data,  acquired  during  July  1984  in  the 
Gulf  of  Genova,  show  that  a  significant  detection  improvement  can  be  achieved  by 
high-resolution  data  processing  for  the  lowest  signal-to-noise  ratio.  This  improve¬ 
ment  ranges  from  5%  to  20%  depending  on  the  allowed  false  alarm  rate  and  on  the 
signal-to-noise  level  considered.  For  the  highest  signal-to-noise  ratio  the  performance 
of  the  two  beamformers  was  comparable. 

These  results  suggest  that  for  acoustic  fields  dominated  by  a  multi-source  noise  struc¬ 
ture,  the  U9e  of  a  high-resolution  beamformer  may  lead  to  an  improvement  of  the 
detection  performance  of  conventional  sonar  systems  This  detection  improvement 
is  due  to  the  ability  of  the  high-resolution  beamformer  to  resolve  the  signal  in  the 
presence  of  such  a  multi-source  structure.  In  this  case  high- resolution  beamforming 
represents  a  potential  candidate  for  low-level  target  detection  with  short  aperture 
arrays.  A  real  life  example  of  such  a  situation  was  shown  in  this  study  for  an  acoustic 
field  dominated  by  ship  traffic  noise. 

However,  further  work  is  needed  to  enlarge  the  statistical  significance  of  the  results 
obtained  here  to  other  databases  acquired  in  various  sites  and  with  different  ship 
distributions.  Possibly,  well-adapted  high-resolution  techniques  can  be  designed 
to  include  a  priori  noise  field  coherence  characteristics  in  order  to  achieve  higher 
detection  performances  in  ship-induced  noise  fields. 
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Appendix  A 
WBS  algorithm  equation  flowchart 

The  following  notation  i«  used:  AU  quantities  in  capital  letters  are  defined  in  decibels. 
St,  is  the  number  of  beams  on  the  conventional  planewave  beamfornter,  and  St,  is 
the  number  of  beams  on  the  high-resolution  held  estimate;  in  practic  St,  -  4St, 
The  quantity  BW„  is  defined  as  the  nth  conventional  besunwidth. 

Step  l  Gel  the  measurements: 

{B„),  n  =  l, . .  . ,  S),. 


Step  2  Get  array  beampattern: 


(H,),  «  =  1 . Sk 

Step  3  Get  first  field  estimate,  k  =  0: 

F,fc  =  min (B„  -  -  10 logfBWJ] ,  i  = 

Step  4  Get  predicted  beamformer  response  to  field  F* 


Bn(F‘)=  10  lcg^io''** 

i=l 


Step  5  Compute  dB  error: 


*  ,r 

DE(fl,F‘)  =  £|«n-.  Sn(F‘)j 

n=  l 


Step  6  Stop  iteration  test  on  dB  error: 

If  DE(B,  Fk)  <  DEo  then  stop,  otherwise  go  to  Step  7. 

Step  7  Update  field  estimate: 

Fk+I  =  F‘ +  C,(B,Fk),  s=l.....Wk, 
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with  the  correction 


Bn  IMPcn> 

("\(  B  t  F  )  =  u>, - — ri - 


where  run  constant  generally  set  to  1.5  or  2  and 


IMP,  ,  =  10 


u,.  -  inin|l,TIMP,/AVlMP,) , 


where 


TIMP,  =  ^1MP„  . 


and  AVIMP,  is  defined  as  TIMP*  for  a  flat  field,  i.e.  F ,  =  constant. 
Step  S  Stop  test  on  the  number  of  iterations: 

If  k  =  Jb0  then  stop,  otherwise  increment  k  -  k  +  1  and  go  to  Step  i. 
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